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^ ' There is some observational evidence for earlier evolution of clusters of galaxies than 

predicted in the standard ACDM model with a Gaussian primordial density fluctuation 
field, and a low value for the mass variance parameter (o-g). Particularly difficult in this 
model is the interpretation of possible excess CMB anisotropy on cluster scales as due to 
CO ' the Sunyaev-Zeldovich (S-Z) effect. We have calculated S-Z power spectra in the standard 

model, and in two alternative models which predict higher cluster abundance - a model 
with non-Gaussian PDF, and an early dark energy model. As anticipated, the levels of 
S-Z power in the latter two models are significantly higher than in the standard model, 
' and in good agreement with current measurements of CMB anisotropy at high multipole 

values. Our results provide a sufficient basis for testing the viability of the three models 
by future high quality measurements of cluster abundance and the anisotropy induced 
by the S-Z effect. 
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1. Introduction 

The power spectrum of the primary CMB anisotropy falls sharply at mutlipoles 
I > 1000, where 'excess' power is primarily induced by the Sunyaev-Zeldovich (S- 
Z) effect. There are some indications that such excess power was measured by the 
CB1 , ACBAR , and BIMA experiments. Attributing the measured power to the 
S-Z effect requires a higher cluster abundance than predicted in the standard ACDM 
cosmological model, particularly so for the low value of the mass variance parameter, 
o"8 = 0.74+q'qq, deduced from the WMAP 3-year dataPl. An unrealistically high value 
°8 ^ 1 would be required for consistency with the current CMB measurements. This 
possible discrepancy enhances interest in alternative models in which higher cluster 
abundances are predicted. 
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Other evidence for earlier and more abu ndan t cluster population comes from 
radio observations of clusters at high redshiftPG3 anc i larger values of the concen- 
tration parameter and Einstein radii than expected in the standard model , which 
imply earlier formation of high mass (M > 10 15 M Q ) clusters. The cluster angular 
two-point correlation function provides yet another measure of the abundance and 
evolved nature of the population. A recent analysis of Spitzer Space Telescope ob- 
servations^ seems to indicate early formation of very massive galaxies, with number 
densities that are considerably higher than predicted in the standard ACDM model, 
and surprisingly high level of clustering. 

The first of two alternative models, in which clusters are expected to form earlier 
than in the standard model, is an isocurvature CDM with scale-dependent non- 
Gaussian, x m distributed primordial density fluctuation 

field™ 

where m is the 

number of CDM fields added in quadratures to yield the x m distribution. With 
increased number of random primordial density fluctuation (PDF) fields, their sum 
approaches a Gaussian distribution (in accord with the central limit theorem); thus, 
the degree of deviation from a normally distributed PDF is the largest for m = 1. 
The evolution of the large scale structure and pri mary CMB anisotropy in the 
X 2 family of models were explored in several studies ^ ! \ showing explicitly that 
primordial overdensities attain larger amplitudes with higher probabilities than in 
a Gaussian field. Correspondingly, cluste r form earlier and are more abundant, 
thereby enhancing levels of S-Z observable J 1 2 | 1 3 [ j n latter two papers the then 
current WMAP 1-year normalization was used, as = 0.9, and it was shown that 
already with this relatively high value it was difficult to reconcile the CMB power 
excess with the inferred cluster population, if the primordial fluctuation field was 
Gaussian. 

We note that the degree of non-Gaussianity in the explored % 2 models is con- 
sistent with limits set by analyses of the WMAP data^. This is largely due to 
the fact that the latter dataset yields information on scales that are much larger 
than those associated with clusters. On these large scales the low overdensities may 
be indistinguishable from a non-Gaussian distribution when the density field is not 
scale invariant. 

Our previous work focused on the predicted levels of S-Z anisotropy and cluster 
number counts in the Xi modeff^^, showing that S-Z power levels in this model 
are appreciably higher than indicated by current measurements. Here we explore 
predicted S-Z power spectra in the x\ model. 

Temporal variation of the dark energy density in early dark energy (EDE) mod- 
els provides an alternative for generating an enhanced cluster abundance at higher 
rcdshifts. In these models the DE is appreciable already at early epochs and at- 
tains the observationally inferred value at present. The evolution of structure in the 
linea r ' regim e and CMB anisotropy in these models have been explored in sev eral 
wor k» Two specific EDE models have recently been investigated in detaiPl, 
resulting in explicit numerical determination of the linear growth factor of density 
perturbations, critical density for spherical collapse, S c , and overdensity at virializa- 
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tion, S v , quantities needed to calculate the cluster mass function. The non- vanishing 
dark energy component at early times drives an early acceleration phase, implying 
a slower evolution of the linear growth factor and reduced values of S c . Thus, for a 
given value of the mass variance normalization, <tm , the corresponding quantity at 
early times should be larger than what is implied in the ACDM model, and since 
the critical overdensity for collapse at a given redshift is linearly extrapolated from 
an earlier time, the slower evolution of the growth factor in EDE models results in 
a lower 5 C as compared with its value in ACDM model. The reduced 5 C and higher 
<jm obviously yield a more abundant cluster population. 

Here we update the predicted S-Z power spectra in the above two alternative 
cosmological models. In Section 2 we briefly describe the variant of the Press & 
Schechter mass function in the x 2_ distributed PDF, and outline the properties of 
the EDE model adopted in our calculations. Results of the calculations of S-Z power 
spectra in the ACDM, EDE, and non-Gaussian models are presented and compared 
in Section 3, followed by a brief discussion in Section 4. 



2. Calculations 

The calculations of S-Z power spectra requires knowledge and modeling of global, 
large scale, and cluster quantities. To do so we adopt the methodology described 
in several papers^212QI. We refer to the ACDM, EDE, and non-Gaussian models 
as models I, II, and III, respectively. Models I and III were described by uJ^l, 
so our brief discussion here will include only the most essential aspects of these 
models. The EDE model we adopt here is characterized by the density parameter of 
early quintessence Q e = 0.03, and the coefficient of the equation of state parameter 
w = -0.9 at z = 0. The effective coefficient as function of redshift i 

J2H 

W{Z) = l + U log(l + z) ' (1) 



where 

-3w 



log(^)+log(^)' 



(2) 



i(M,z) = -^) T7 ^—dM, (3) 



and Q m is the matter density parameter. 

The basic quantity in our calculations of S-Z power spectra is the Press & 
Schechter (1974) mass function, 

Pb da 
' Ma\ l {z)~d~M 

where p = 8 c (z) / <jm(z) is the ratio of the critical overdensity for collapse to the 
mass variance ctm at redshift z, and pt is the background density at z = 0. For a 
Gaussian PDF field assumed in models I and II, we have 

^) = \/fe-^ 2 / 2 >. (4) 
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The corresponding expression for a PDF which obeys x| statistics of model III is 
particularly simple 

F( M ) = e-( 1+ ^. (5) 

The mass function is normalized such that all the mass is included in halos, a 
normalization that in the original Press & Schechter mass function was affected 
by including a ('fudge') factor of 2. The functional form of the mass function is 
different in model III from that in models I and II, and the redshift dependence of 
the critical density for collapse and the linear growth factor both differ in model II 
than the corresponding quantities in the oth er tw o models. Explicit expressions for 
these quantities are given in our recent papei^. 

The mass variance <tm was calculated with a top - ha t window function and CDM 
transfer functions taken from Bardeen et al. (1986)^ - adiabatic transfer function 
for models I and II, and isocurvature transfer function for model III. The shape 
of the CDM transfer function in the EDE model is slightly different than in the 
standard model; this difference is ignored here. For the isocurvature model this 
function (of the wavenumber k) is 

1 -5/4 



T{k) lsoc = {5.6qY 



1 + Mr + (5 6q)^ 

1 + 215g + (16g) 2 (l + O.Bq)- 1 1 q> 



(6) 



where q = k/ \fl m h 2 Mpc" 1 ), with h = i? /(100 km s" 1 Mpc" 1 ). Values of the global 
parameters were taken to be those deduced from the WMAP 3-year data, f^A = 0.76 
(which, for the case w ^ — 1 is usually written as J1q), £! m = 0.24, h = 0.73, and 
erg = 0.74^Q Qg. The spectral index of the PDF spectrum is n = 1 in models I and 
II, and n = — 1.8 in model III. In the EDE model the differential equations for the 
evolution of the linear growth factor, 5 C , and 5 V , are different than in the other two 
models. Full description of the of these equations and their numerical solutions can 
be found in our recent papeJ^. 

A meaningful comparison of the predicted S-Z spectra of the three models con- 
sidered here requires proper and self-consistent normalization of the respective mass 
functions. This is accomplished by requiring that the cumulative cluster density at 
z = is in agreement with that calculated in the standard ACDM model. 

In addition to the global and large-scale parameters, the calculation of S-Z power 
spectra necessitates full description of the properties of IC gas. To do so clusters 
are assumed to be in hydrostatic equilibrium with [3 profiles for the total mass and 
gas density profiles, and with a gas mass fraction of 0.1. The gas temperature is 
determined from the equation of hydrostatic equilibrium. 



3. Sample Power Spectra 

The behavior of the cumulative mass function in the three models at z = 3 and 
z = 0.01 is shown in Fig.[T]for the mass range IO^Mq/i" 1 < M < IO^Mq/i^ 1 . As 
is clear, the models are correctly normalized to yield the same cumulative cluster 
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Fig. 1. Cumulative mass functions in the three models. The upper and lower curves for each 
model correspond to redshifts z = 0.01 and z = 3, respectively. 



density at low redshifts. Abundances of high-mass clusters are indeed higher in 
models II and III than in model I. As anticipated, the relative abundances with 
respect to those in ACDM increase with redshift. It is also apparent that at early 
times the excess of high overdensity fluctuations in the xi m odel has a stronger 
impact on the abundance of high mass clusters than the slower evolution of the 
linear growth factor and lower values of S c in model II, whereas their respective 
effects at present are about the same in these two models. 

The enhanced cluster abundances in models II and III are directly manifested 
in higher S-Z power levels than in the standard model, as is immediately evident 
in Fig. [2] which shows the spectra at 31 GHz together with CBI, ACBAR, and 
BIMA measurements at this frequency. In addition, the (broad) peak power in the 
non-Gaussian model is reached at multipoles, I ~ 7000 — 8000 as compared with 
I ~ 4000 — 5000 in model I and II, a direct consequence of the higher abundances 
of distant clusters (with smaller apparent sizes) in model III. 

The predictions of both the EDE and non-Gaussian models are more consistent 
with current measurements of the CMB power spectrum. The shape of their S-Z 
spectra are virtually identical for I > 10 4 , with power levels somewhat lower at 
lower £ in the former model than in the non-Gaussian model. 
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Fig. 2. S-Z angular power spectrum at v = 31GHz in models I (continuous), II (dashed), III 
(dash-dotted), and IV (thick dash-dotted). The shaded areas correspond to the WMAP reported 
1-(T errors in erg. Also shown is the power excess reported by the BIMA (diamonds), CBI (crosses) 
and ACBAR (x-symbols) experiments. Note that the BIMA indicated power at I = 8748 is an 
upper limit. 

4. DISCUSSION 

The work reviewed here has been motivated by initial indications that there might be 
a significant discrepancy between current measurements of levels of CMB anisotropy 
on scales t > 2000, and predicted levels of S-Z power in the standard ACDM with 
gaussian PDF field. The discrepancy stems mainly from the fact that high-mass 
clusters, the largest contributors to S-Z power whose density decreases sharply with 
decreasing erg , are not sufficiently abundant if this important parameter is as low as 
deduced from the 3-year WMAP data. A higher cluster abundance is expected when 
the PDF field has the form predicted in the non-gaussian xl model. This is a result 
of higher probabilities for overdense regions at high z, leading to earlier collapse 
of proto-cluster halos. In the early quintessence model earlier cluster collapse is a 
manifestation of higher linear growth factor, and lower value of the critical density 
for collapse at high z. Accordingly, levels of S-Z power are higher in these two 
alternative models. Clearly, these two non-standard models are by no means the only 
viable alternatives to the standard model. As we noted, the non-Gaussian model 
considered here is just one of the x m family, but with decreasing non-Gaussianity 
with increasing m. Also, other early quintessence models with higher EDE densities 
result in slower evolution of the linear growth factor, and reduced value of S c , thereby 
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leading to higher levels of S-Z power. 

While levels of S-Z power spectra predicted in the standard model span an ap- 
preciable range, reflecting also uncertainties in the evolution of internal properties 
of clusters - such as IC gas density and temperature - their maximal values are still 
well below current observational results. This concl usion i s based on extensive inves- 
tigations of the cluster temperature-mass 

relationP323l 

IC gas models, including 

non-isothermal polytropic temperature profiles™^ 1 , and evolution of the gas mass 
fraction. 

Our purpose here has been to show that the two alternative models explored 
can produce S-Z power levels that are substantially higher than in the standard 
model. It is too early to actually fit the predictions to the preliminary high £ re- 
sults. This will have to be done in conjunction with other cluster observables, such 
as cluster (S-Z) number counts, and the two-point correlation function (which were 
considered in our previous worlffSl). Other cluster measures that are very much 
affected in these models are formation times and concentration parameters (whose 
observational manifestations include, e.g., mass profiles and Einstein ring sizes). 
Nonetheless, it is quite clear from our results that the xi model, and EDE models 
with larger quintessence densities at early times than in the specific model adopted 
here, predict significantly higher levels of S-Z power than indicated by current ob- 
servational results, and therefore seem to be non-viable. 

We should know soon whether the apparent discrepancy is real, when results of 
more extensive and precise measurements of the high £ power will become known. 
But irrespective of these upcoming CMB and S-Z measurements, there seem to be 
other observational indications that cluster formed earlier and are more abundant 
than predicted in standard ACDM. The non-gaussian and EDE models explored 
here seem to be viable alternatives, should the discrepancy persist. 
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